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ABSTRACT 

This study establishes the existence of damping structu- 
ral properties in the San Andreas Fault zone through an input 
model of surface loading. Any changes in the density, mass 
or bulk elastic properties of a fault zone will not only pro- 
duce a change in the gravity field of the earth, but also a 
change in the earth's response to the tidal field generated 
by the moon. Knowledge of these changes is important to en- 
gineering constructions, pre-programmed defense systems, etc. 
To adequately measure these responses, sites on opposite 
sides of the San Andreas Fault zone were selected and tidal 
gravity was continuously measured for finite periods of time. 
The Pacific Ocean and its very low frequency lunar portion of 
the ocean load spectrum was used as a driving function to 
compare observed data with theoretically generated data for 
the same sites. A least squares criteria was used to fit the 
observed and theoretical data to the nine most important 
spectral lines and to obtain the phase and amplitude for each 
line at both sites. Utilizing the most stable lines to com- 
pare values, a significant drop (about 10$) in both phase 
and amplitude occurred between the two sites. The anomalous 
decrease in both areas of analysis discloses a structural 
damping in the San Andreas Fault zone. 
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I. INTRODUCTION 



A. STATEMENT OF THE PROBLEM 

The combined gravitational fields of the moon and the 
sun acting on the solid earth are sufficient to displace its 
surface on the order of 20 cm. semi-daily [ 1 ] . Unfortunate- 
ly, this direct forcing function is so broad in its areal 
extent (i.e., several thousand kilometers) that it cannot 
reasonably be used to study responses of small crustal 
structures (such as fault zones) [2] which are on the order 
of a few kilometers in width and depth. However, the response 
of the ocean to the lunar-solar gravitational field (on the 
order of meters) constitutes a sufficiently energetic forc- 
ing function such that under certain conditions dynamic 
properties of crustal structures can be studied. 

The conditions under which the ocean can effectively be 
used as a forcing function to study the properties of a 
shallow structure are defined as: 

1. The distance between the edge of the ocean load and 
the fault structure should not exceed the smallest dimension 
of the fault zone by one order of magnitude [3]; 

2. Noise response in the earth tide records obtained 
during this experiment must be kept to a minimum. Minimiz- 
ing noise response is achieved by using parts of the ocean 
tidal spectrum which are both stationary and relatively un- 
affected by meteorological influences. 



The two lines chosen for amplitude and phase analysis in 
this study are the principal lunar diurnal 0-^ and the princi- 
pal lunar semi-diurnal lines which are the most energetic 
in the theoretical tidal spectrum and are the most removed 
from possible contamination due to meteorologically induced 
noises. 

Earthquake prediction research explores the causes of 
earthquakes from two points of view: (1) failure brought 

about by an external increase of shear stress which exceeds 
the shear modulus of the fault; (2) failure brought about 
when the internal properties of the fault change (coefficient 
of shear) until they are reduced below that of the external 
shear being applied. 

In this study responses to tidal gravity were measured 
at two points on either side of the San Andreas Fault. Tidal 
gravity is a vertical measure of the acceleration due to 

gravity. The value of gravity is on the order of 10 p cm/sec/~ 

3 —6 

sec or ICr gal. Resolution on the order of 10 gal is 

achieved with the instruments used (described below) giving 

Q 

results accurate to one part in a 10 . The distance between 
the two stations (about 5 kilometers) is small compared to 
the distance to the edge of the ocean load (Pacific Ocean, 
about 30 km). The amplitude and phase of the two measured 
lunar tidal constituents (0-^ and N^) were compared with the 
amplitude and phase of the same two lines derived from an 
analysis of the theoretical tidal gravity calculated at each 
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of the two sites for identical periods of time. The differ- 
ence in the observed phases between the sites was compared 
with the difference in theoretical phases between the sites 
and a 10 $ change in the means was observed. 

This result suggests that the principal source for the 
phase difference between the observed and theoretical values 
at both sites is due to the response of these sites to ocean 
tidal loading at the same frequency. However, a 10$ devia- 
tion of the difference between these two sites from the mean 
is significantly large and is not due to a simple shift of a 
few kilometers further away from the coastline. Instead, the 
anomalously large difference in the comparison is explained 
by the fact that the stations are not just separated by a 
few kilometers but that they lie on opposite sides of a weak- 
ness in the earth's crust. The absolute difference between 
the most inland station and the other is a lag of about 
three degrees at the Mg frequency. This is interprted pri- 
marily as a viscoelastic response of the fault zone to ocean 
tidal loading. The inland edge of the San Andreas Fault 
reaches it peak amplitude of ocean tidal response approxi- 
mately 1/10 hour after the peak response of the ocean side 
of the fault. 

While this study did not attempt to directly calculate 
any geological strength parameters, it did suggest the 
existence of a crustal weakness with viscoelastic properties 
that strongly dampen the Mg line of the ocean tidal loading. 
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B. BACKGROUND 



There exists in any multi-body system a measurable gravi- 
tational potential due to the attacting fields of the many 
masses (or bodies) present. Since all bodies have motion 
relative to all others in a system, the potential fields 
fluctuate with time. On the earth, only the sun and the moon 
contain sufficiently large quantities of mass at a close 
enough distance to significantly affect tidal variations or 
time fluctuations in the potential field. 

The effects from the sun and the moon generate earth 
tides in a similar manner to their production of ocean tides. 
Unfortunately , the study of earth tides is complicated by 
the presence of oceans which also cause changes in earth 
tides. Although quite small, ocean loading effects on the 
earth tides have been detected in measurements as far inland 
as mid-continent in the United States [4]« The mass of the 
oceans complicates measurements of the earth tides in the 
same ways as external masses (the sun and moon) plus (1) 
the loading effect due to water (additional mass) and (2) 
deformation of the earth due to the redistribution of ocean 
waters as 'the water more freely fluctuates. 

The combination of effects of the sun, moon and oceans 
over time results in a continuous change in the earth's 
gravity. The gravitational force at any given point in time 
and space can be broken into two components, one normal to 
the plane tangent to the earth's surface and one in the 
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plane tangent to the surface. This study will concentrate 
on the former component called tidal gravity. 

Past geological and seismological studies have examined 
the bulk elastic properties of the whole earth and the dis- 
tribution of these properties in the core and the mantle. 

Few studies have narrowed their objective to examining the 
bulk elastic properties of only the crust of the earth and 
in a relatively small area. Seismic studies show that the 
San Andreas Fault zone has brittle properties to a depth of 
only ten kilometers, which is neither the bottom of the 
fault nor the bottom of the continental plate [5]. It is 
approximately one-third of the depth of the plate. Authors 
of these studies quite often generalize their results to 
large area averages and tend to overshadow the relatively 
small area of a shallow fault zone in the final analysis [2]. 
The seismic studies and their related high frequencies (from 
earthquakes) also do not allow determination of anelastic 
responses, but very low frequency studies such as those re- 
sulting from ocean loading should disclose any anelastic 
properties present. 

This study is a portion of a larger project, (headed by 
Dr. M. D. Wood) at the National Center for Earthquake Re- 
search (NCER), to establish that the response of the earth 
to tidal loading is significantly perturbed by the San 
Andreas Fault and also to determine through spatial and 
temporal means the bulk properties of the fault zone. 



12 



There are two possible sources of error in analyzing 
these data that should be noted. First, due to many con- 
straints, the three experiments in this study were performed 
sequentially instead of simultaneously. Second, the San 
Francisco Bay, due to its size and its proximity to the sites, 
will affect the values measured in the experiment. The ef- 
fects of these errors are assumed to be relatively small, 
but their exact effects are undetermined. 
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II. DATA ACQUISITION AND PROCESSING 



A. LOCATION SELECTION 

The locations selected for the two main sites for this 
analysis were obtained from a threefold criterion. First, 
it was desired to ensure that the two sites were on opposite 
sides of the San Andreas Fault zone and nearly equidistant 
from the fault region. This was necessary in order to ob- 
tain spatial analysis of the tidal wave as it passed from one 
site to the other by traversing through the fault region. 
Second, the sites were desired to be on a line nearly per- 
pendicular to the coastline. Third was the availability of 
suitable sites. It was imperative that the physical location 
of each instrument be based on a very sturdy foundation, 
ideally on deep lying rock, to minimize noise. To further 
eliminate noises (in the form of other energies such as wind, 
temperature fluctuation, surface vibrations, etc.) it was 
desirable to place each recording instrument below the sur- 
face of the earth. Time and financial restrictions forced a 
compromise of this desire to a level, surface (or very 
shallow) base that was protected from as many of the other 
noise-producing sources as possible. In addition to these 
restrictions, the security of the expensive recording instru- 
ments became a major concern. The instruments could not be 
left in the open where vandals or sightseers might tamper 
with them. It was also desired to have daily time marks 
entered on each record, generating the requirement for 
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someone to be present at the site each day to perform this 
task. Appropriate sites with the necessary attendance re- 
quirements we re located [Figure 1]. One recording instru- 
ment was in a basement vault on a base on concrete poured on 
solid rock (Block experiment); another was in a separate 
garage on a concrete base poured on rock (Rose experiment). 
Both properties were private and actual instrument locations 
were secure. The third site was located in the NCER Build- 
ing in Menlo Park on a concrete floor (NCER experiment). 

B. EARTH TIDE RECORDINGS 

All the earth tide recordings were taken from LaCoste- 
Romberg Earth Tide Gravity Meter type instruments [l] on 
paper strip charts [Figure 2]. These were recorded from 
NCER owned instruments run by 120v AC with l2v automobile 
batteries as a back-up. 

To allow future adjustments for inaccuracies in chart 
speed, stretching or shrinking, etc., daily time marks were 
mechanically entered on each chart and labeled. Chart speed 
for each experiment was two inches per hour with a displace- 
ment of one inch equal to 30 y^gal. Each experiment was 
attempted to run at least thirty days to give a large enough 
sample for adequate analysis [6]. 

C . DATA PROCESSING 

The strip charts were digitized onto magnetic tape by 
use of a Calma 302 digitizer. The cross-hair cursor moved 
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Figure 1. i Test Site Locations. 
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Figure 2. Raw Data of Earth Tide from Rose Experiment. 



vertically or horizontally as desired across the strip 
chart taking discrete samples every .01 inches and recording 
the readings. A special indication was noted at each daily 
time mark entered on the chart by monitoring personnel, and 
the time and date were recorded by the person performing the 
digitization. The digital tape data was later corrected 
for chart stretching, shrinkage, etc., by proportionally 
distributing the actual recorded values between time marks 
to the correct number of .01 inch increments (.01 inch 
corresponds to 1$ seconds of time) required between the daily 
time marks. This localized physical errors between time 
marks (usually every 24 hours) and minimized accumulation 
of error. 

The digitized data consisted of readings based on a 
continuous time record. These data were then processed with 
a computer program (designed by M. D. Wood and P. A. Marshall 
[7] from the National Center for Earthquake Research (NCER) 
and modified by M. D. Wood and R. L. Keller for use in this 
study) which decimated the sample interval from .01 inches 
to hourly values taken on the hour. 

Sampling theory requires a sample interval at least twice 
as often as the highest expected frequency. Implicit in this 
is the assumption that the recorded data be of a reasonable 
duration of time. A reasonable time duration is propor- 
tional to the reciprocal of the smallest difference between 
any two frequencies expected in the data [6]. Since only 
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nine lines of the BO-plus theoretical lines are used, it can 
be calculated that at least one month of data is necessary 
for a 1 io resolution of tidal data [6]. The smallest fre- 
quency line used in this study is the with a frequency 
of less than 4 cycles/day, requiring at least 8 samples per 
day. By oversampling, the resolution remains well within ifo. 

III. LEAST SQUARES ANALYSIS 
A. AMPLITUDE ANALYSIS 

The observed tidal gravity records were fit with a least 
squares computer program to the nine theoretical lines of the 
earth tide spectrum (explained in Appendix A). The computer 
program (designed by M. D. Wood [6] and modified by M. D. 

Wood and N. E. King of NCER and R. L. Keller) performs a re- 
gression of the data to fit both an n— polynomial and a 
sine wave function simultaneously. The results give the 
amplitudes and phases of the nine lines mentioned above. 

It is a well established fact that the principal lunar 
semi-diurnal M^ and the principal lunar diurnal 0-^ lines of 
the earth tide frequency spectrum are the most invariant of 
all the tidal lines [l], [6], [8]. Of these two, the M 2 
line has proven to be less variant from its theoretical value 
than the other lines in previous studies conducted in this 
region [6], [7], [8] and was thus used as a normalizing fac- 
tor for all of the spectral lines. The normalized amplitudes 
were computed using the following formula: 
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A j/ A M2 



where is the (dimensionless) relative amplitude, is 
the amplitude of the constituent line being normalized, and 
A^ is the amplitude of the M2 line. This normalizing pro- 
cess results in dimensionless values, facilitating compari- 
son between observed and theoretical results. 

The normalizing of the amplitudes of the selected 
lines in the earth tide spectrum was conducted for both the 
theoretical data [Table III] and for the observed data [Table 
IV]. Visual displays of these data are given in Figure 7 
and Figure S. 

Taking the expected amplitude (theoretical) of each 
constituent line for each of the sites and subtracting the 
corresponding observed amplitude, the amplitude differences 
are obtained [Table I]. Differences can be expected due to 
forcing nine lines of a frequency spectrum to absorb what is 
acknowledged to consist of at least $9 lines and also due to 
other non-tidal factors such as seasonal noise, extreme wea- 
ther conditions, etc. These differences should remain 
about constant between sites (resulting in low values in the 
last column) and would indicate a high stability of a con- 
stituent line. The last column of Table I shows a standard 
deviation of about 5$ for the 0-^ line from its average theo- 
retical mean. This indicates that the 0-^ line could be pre- 
dicted to have a measured (observed) value with a standard 
deviation equal to about 5$ of the theoretical value of the 
amplitude. 
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The normalizing of the amplitudes of the selected 
lines in the earth tide spectrum was conducted for both the 
theoretical data [Table II] and for the observed data [Table 
III]. Visual displays of these data are given in Figure 3 
and Figure 4. These tables and figures are presented here 
to facilitate comparison of the results obtained in this 
study to those of other completed studies [6], [7], [8] 
making use of similar types of displays. 

A relative gravimetric factor G^f can be computed by 
the following formula: 

( experimental )/A^ (theoretical) 
r* 

where the A^ are described above. Using the two most stable 
lines to establish a comparative base, we compute: 

G r (0-^/M2) = A. r )( experimental )/A r ( O-^/M^ ) (theoretical ) 

Theoretically, this ratio should equal to 1.00, but with 
other dynamic considerations a 1$ variation [9] is accepted. 
Ten world-wide stations studied by Harrison and others [10] 
resulted in an average value of : 

G^Ol/M^ = 0.939+0.035. 

Four stations in the San Francisco Bay area analyzed by Wood 
[6] gave an average value of: 

G r (0 ] _A 2 ) = 0.938+0.026. 
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The average value obtained by the three stations in this 
study gave: 

G r (0 1 /M 2 ) = 1.051+0.061. 

These values should not be different, regardless of loca 
tion, yet this study and that by Wood indicate values about 
4$ to 5 i° outside of the accepted range. Both of these 
studies were conducted on the West Coast near the San 
Francisco Bay Region and could indicate a local disparity 
from areas with thicker crustal depths (The California coast 
al regions are known to have a much thinner crustal depth 
than most other continental areas). However, no conclusion 
is drawn from this study other than to mention that this dis 
parity should not exist. 

B. PHASE ANALYSIS 

The observed tidal gravity records were simultaneously 
fit in the least squares program (described in the previous 
section) to the nine theoretical waves of the earth tide 
spectrum to estimate the phase of each constituent line. 

The values for both theoretical and observed phases for the 
M 2 lines are given in Table IV. All charts were maintained 
in reference to Greenwich Time (GMT) except for the Rose 
site which was marked in local daylight time. A seven hour 
difference existed between these base reference times and 
the phase was appropriately advanced for this discrepancy. 
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The results shov/ a phase difference between the theoreti- 
cal and observed line at both the Rose site and the Block 
site. 






TABLE I. RELATIVE EARTH TIDE FACTORS - COMPARATIVE RESULTS 
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TABLE II. RELATIVE EARTH TIDE AMPLITUDE FACTORS - THEORETICAL 
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TABLE III. RELATIVE EARTH TIDE AMPLITUDE FACTORS - OBSERVED 
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RELATIVE AMPLITUDE 



R = ROSE 



B = BLOCK 



N = NCER 




Figure 3- Composite Theoretical Amplitude Histogram 
Normalized by Mg Component 
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RELATIVE AMPLITUDE 



i 



R = ROSE 
B = BLOCK 




Figure l+. Composite Observed Amplitude Histogram, 
Normalized by Mp Component. 
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TABLE IV. RELATIVE EARTH TIDE PHASE COMPARISONS 
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FORTNIGHTLY WAVE 
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IV. RESULTS 



The average value of G^O-^/Mg) was much higher than those 
produced by the experiments done by Harrison and others [10] 
and Wood [7]* As noted above, no explanation for this has 
to date been found. The relative gravimetric factor showed 
a downward change between the two sites. It decreased by 
more than 0.1, or about a 10$ drop in value. 

In the examination of the two most stable lines in the 
earth tide spectrum, the value of the 0-^/Mg (the normalized 
amplitude of the 0^ constituent line) decreased from the Rose 
to the Block station by a value of .OUBB which is about a 
7 4$ decrease between the two sites. Theoretically, the nor- 
malized 0-^ line amplitude should not show much change in 
value, specifically across a distance of 5 km. From the 
values in Table II, it can be shown that about a 2$ increase 
should have occurred. 

There also exists a phase difference between the theore- 
tical and observed data at each site. The observed and 
theoretical lines should all be in phase with each other. 
Allowing for small variations due to noise and tectonic in- 
terferences prior to the arrival of the Mg line at the Rose 
site (again the most stable line is selected for comparison), 
there may exist a phase difference between theoretical and 
observed lines. Between the Rose and Block sites, however, 
this difference should remain constant. In fact, it does 
not; it changes by more than 2-g- degrees. This change is a 
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decrease in phase difference, or a lag, of about 9ir i° between 
the two sites. Physically interpreted, this is a time delay 
of about 1/10 hour for peak amplitude of the wave to 
travel from the Rose to the Block site. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



With the evidence advanced in this study, it appears 
that some disturbing factor has caused a large percentage 
decrease in tidal amplitude and phase between the Rose and 
the Block sites. These percentage changes are much too great 
to be explained by the 5 km. distance between the sites. The 
fact that the San Andreas Fault is located between these two 
sites leads one to suspect that drastic changes in earth 
tides between the two sites might be due to the fault region 
and its different structural properties from surrounding area 
properties . 

Assuming that the fault region is responsible for these 
tidal changes, it is possible to calculate the approximate 
strength parameters of the geological structure of the fault 
region itself through a more detailed geological study of 
sites across the fault. 

Since this experiment was limited in time, money, equip- 
ment, available sites, etc., it cannot provide absolute 
proof that the fault region is responsible for the observed 
changes. The study provides evidence that the fault might be 
reason for the observed changes. It is recommended that 
further investigations of a similar nature be conducted at 
sites void of an inland bay to verify the res\xlts of this 
study. These future studies could also be set up to simul- 
taneously examine the specific geological properties exhibited 
by the fault region in each area of study. 
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APPENDIX A 



THEORY OF EARTH TIDES 

A study of the earth tide curve shown in Figure 5 reveals 
the presence of at least two frequency components of the 
earth tide spectrum. Identified in this figure are the diur- 
nal frequency with a period of approximately 28 days and the 
semi-diurnal frequency with an approximate 14 day period. A 
Fourier analysis of a record such as this reveals that there 
exist many constituent waves that combine to give us the wave 
form. Analysis of longer records would reveal semi-annual, 
annual, and even longer period waves which correspond to mo- 
tions of the sun and moon in relation to the earth. 

When the complex relationships between the moon (and sun) 
of a rigid earth are considered, the constituent waves (lines 
in the spectrum) derived via the Fourier analysis become 
meaningful. Let us examine a fictitious model [6], [S] in 
which the moon is orbiting about the earth in the plane of 
the earth's equator at a constant distance with a constant 
velocity equal to the mean distance and velocity of the actu- 
al moon. The earth is also revolving on its axis once every 
24 hours and this motion plus the moon's gravitational at- 
traction on the earth's surface produces a tide with a per- 
iod equal to a lunar half-day and is described by the M 2 
principal lunar semi-diurnal line, the most energetic of the 
semi-diurnal lines. In similar fashion, the S 2 principal 



34 



solar semi-diurnal line is produced by a fictitious sun in a 
circular orbit at a constant velocity and distance. 

If the moon's declination in its orbit (above or below 
the plane of the earth's equator) is now added to our simple 
model, the distance between the moon and an observation point 
on the earth varies with time, but in a different manner than 
the simpler model. In the simple model, the distance varia- 
tion with time was due only to the observer's position and 
the lunar position in the equatorial plane. Now a declina- 
tion distance away from the equatorial plane also varies with 
time as well as the consideration for the earlier model, add- 
ing another variable to be considered in the distance varia- 
tion between the moon and an observation point on earth. 

For observation points on the equator, the effects of the 
lunar orbit in the northern and southern hemisphere are equal 
in magnitude. For non-equatorial observers, the amplitude 
of -the two diurnal maxima would vary in time with the chang- 
ing distance from the observation point to the moon. This 
is referred to as the diurnal inequality. The cyclical mo- 
tion of the moon's declination causes a fortnightly modula- 
tion of the earth tides. 

Adding the ellipticity of the moon's orbit to our model 
now leads to the cyclical changes in the distance from the 
moon to the earth and causes modulation of both the amplitudes 
and phases of all the above effects. 
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Associating a theoretical wave with each of the motions 
of the moon does not complete an explanation of the earth 
tide spectrum. Secondary waves are apparent due to the beat 
frequencies of the principal tide waves. These beat fre- 
quencies are obtained from the sum and difference frequencies 
of the two principal waves. When the difference in the lunar 
and solar principal semi-diurnal waves (M2 and S2) is taken, 
a fortnightly modulation occurs at the beat frequency. The 
lunar elliptic semi-diurnal wave, N2 > and the M2 wave com- 
bine for a beat frequency to form a modulation with a monthly 
period that is associated with the eccentricity of the lunar 
orbit. The luni-solar semi-diurnal declinational wave is 
composed of lunar and solar components which are analytically 
inseparable, as is the K-^ luni-solar diurnal declination wave. 

Further beat frequencies exist which yield waves derived 
in a similar fashion. Table I lists the principal theoreti- 
cal- tidal waves of the earth tide spectrum which are of 
importance to this study. The amplitude of each of these 
waves is dependent on latitude and varies with time. 
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APPENDIX B 



THEORETICAL TIDE GENERATION 



The gravitational potential at some point P on the sur- 
face of the earth [Figure 6] due to another mass (moon or 

sun) can be expressed as [6], [11], [12]: 

Oo 



W = 



0 £2 / \ oo 

!E(s) n . P n (z) = 

^ n=o n=o 



W 



n 



where G is a gravitational constant 
M is the external mass 

r is the radius of the earth 

R is the absolute distance from the earth 
to the external mass 
P n (z) is the Legendre polynomial 

z is the angular position of the external 
mass relative to the zenith of point P 
P is a point on the earth's surface. 



Expanding the terms in the summation above necessitates the 



calculation of the Legendre polynomials. The first few terms 
of the Legendre polynomials are [12]: 

p o< 2) = 1 

P 1 (z) = z 

P 2 (z) = -H3z 2 -1) 

P 3 (z) = K5z 3 -3 ) 

P 4 (z) = l/B(35z 4 -30z 2 +3) 



Examining these results it is seen that the first term, 

W =GM(1/R) , is constant and has no effect on changing the gravi- 
tational potential. The second term, W-^--GM(r/R )cos(z), uni- 
formly affects all points on the surface of the earth and has 
no effect on changing the potential. The third term (and all 
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remaining terms) affect the tidal potential and are all 
candidates for significant effects in the summation. The 
r/R term in the summation is approximately l/60 for the moon 
and 1/2400 for the sun and after n=3 accounts for less than 
2 $ of the previous term. Thus, the gravitational potential 

can be reduced to the simple expression: 

2 2 
X* I* 

W 2 = GMJJ P 2 (z) = -|(3cos 2 (z)-1). 

To express z in terms of fundamental angular parameters 
of position, spherical trigonometry is applied where: 



6 = declination of external mass M 
t = hour angle of point P 

= zenith angle of external mass at point P 
f = latitude of P 
X = longitude of P 
z = polar angle of external mass M 
Using the trigonometric relation: 

cos(cc)= cos(e)cos(z) + sin(e)sin( z)cos(t ) 
and cos(e)= sin(? ), sin(©) = cos (?) 

cos(z )= sin( 6 ), sin(z ) = cos (5 ) 

the gravitational potential becomes 

2 9 
W 2 =GM^j •§■ (3[sin( ? ) sin( 6 ) + cos(?)cos(6 )cos(t)] -1) 
R 



and simplifies to 

W 2 = 3K(sin 2 (?)-l/3)(sin 2 (6 )-l/3) 

+ K(sin(2'f)sin(26)cos(t ) ) 

+ K(cos 2 (£)cos 2 (<$ )cos(2t) ) 



where K = 3/4GM-^ 

R ' 5 



38 






I 















This equation is shown in Figure 7. 

Using the moon as an example, the first term in this ex- 
pression contains the variable <5 ( is constant for a given 
point P) , namely sin (£) , which is an even function varying 
as the declination of the mass M traverses from its maximum 
positive value to its maximum negative value and back each 
month. The functional value stays positive (since sin (5) 
is an even function) and the frequency of the first term be- 
comes two cycles per month and is called the fortnightly wave. 
The second term varies with S and t, but S completes a cycle 
in one month while t has a period of about one day. There- 
fore, this term has a frequency of one cycle per lunar day 
and is called the diurnal wave. Similarly, the third term, 
with a frequency of two cycles per lunar day, is called the 
semi-diurnal wave. 

A similar situation occurs with the sun as the external 
mass, and the sum of the results from the moon and the sun 
yields the terms of the tidal potential. The remainder of 
the tidal potential spectrum consists of the higher ordered 
terms and combined lunar-solar terms neglected in this analy- 
sis. For the purposes of this study the nine waves described 
in Table V are sufficient since they contain more than 90$ 
of the effects of the 39 theoretical lines. 

In the generation of the theoretical earth tides for 
this study, J. C. Harrison's program [13] was used. In the 
development of his model for generating earth tides, Harrison 
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uses W2 for the moon (as described above) and W2 for the sun 
and also includes for the moon to achieve greater preci- 
sion in his results. 
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POLE 




Figure 6. Geometric Relation between Station Coordinates and 
Tide Producing Body Angular Coordinates. Where ct 
is the zenith angle of the tide producing body at 
the station location P. P is at colatitude S and 
longitude A . Z is the polar angle and 1? the 
terrestrial east longitude of M. (Adopted from 
Munk and Cartwright ) 
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Figure 7 . AMPLITUDE SPECTRA AND GEOMETRICAL • CHARACTERISTICS OF EARTH TIDES 

(Adopted from Melchior [ 1 ] ) 



TABLE V. TIDAL CONSTITUENTS USED IN ANALYSIS 
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